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Abstract 

We show that the one-dimensional lattice model proposed by Lipatov to describe the 
high energy scattering of hadrons in multicolor QCD is completely integrable. We identify 
this model as the XXX Heisenberg chain of noncompact spin s = and find the conservation 
laws of the model. A generalized Bethe ansatz is developed for the diagonalization of the 
hamiltonian and for the calculation of hadron-hadron scattering amplitude. 
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1. Introduction 



Recently it was suggested by Lipatov |T| that the asymptotic behavior of the hadron-hadron scattering ampli- 
tudes in QCD in the limit of large invariant energy s and fixed transferred momentum t can be described by 
means of the quantum inverse scattering method [| . In the generalized leading logarithmic approximation 
O H' El the scattering amplitude is given by the contributions of the Feynman diagrams having a conserved 
number n of reggeized gluons in the t— channel |Q ||. For fixed n the scattering amplitude can be decom- 
posed into t— channel partial waves f u , describing n gluon to n gluon scattering with total angular momentum 
j = 1 + u). The partial wave satisfies the Bcthc-Salpctcr equation whose solutions can be expressed in terms of 
the wave functions \ 01 the compound states of n gluons as follows jl], [| : 



U{bi},{b^})=J2 J dH ° Xa({£};te($};S)) 



Here, Xa{{bi}; bo) is the wave function of the compound state of n reggeized gluons with transverse coordinates 
{bi} = b±, . . ., b n and it is parameterized by quantum numbers a and by two-dimensional vector bo which 
represents the center of mass coordinate of the states. In the special case n = 2, the function x(bi, b^\ bo) is the 
wave function of the BFKL pomeron jlO[ [l2|. Let us define for all impact vectors bi = (xi,yi), i = 1, . . . , n, the 
following complex coordinates: Zi — Xi + iyi and Zi — Xi — iyi and similarly for the vectors b[. 

Then, in the large N limit, for an arbitrary number n of gluons, the wave functions </? are the eigenstates of 
the holomorphic and antiholomorphic hamiltonians Q|, 11, ^) 

HnX(i z i, Zi}; z , z ) = e n x({zi, Zi}; z , z ), Hnxifa, z^; z , z ) = £ n Xu>({ z h 2 i}; z , z~o), (1) 

with eigenvalues related to complex angular momentum via u> = —j^N(e n + e n ). The hamiltonians H n and 
H n describe nearest neighbors interaction of n particles 

n n 

Hn = / ] Hk,k+i, H n = ^^Hk,k+i (2) 

k=l fc=l 

with periodic boundary conditions H n . n +\ = H n .i, H n .n+i = H n ,i and the two-particle hamiltonians are giving 
by the equivalent representations 

Hj, k = Pr 1 log(z 3 -z k )P J +p- 1 log{z J -z k )P k +log{P J P k ) + 2 lE 

= 2 \og( Zj - z k ) + ( Zj - z k ) logiPjPk^Zj - z k )~ l + 2 lE (3) 
where Pj = ijjj- and is the Euler constant. The same operator can also be represented as jllj 

Thus, in the large N limit the original problem of the calculation of partial waves is reduced to the diagonalization 
of two hamiltonians corresponding to spin chain models defined on a one-dimensional lattice with the number 
of sites equal to the number of reggeized gluons in the t— channel. However, there are additional constraints 
on the eigenstates of the QCD hamiltonians imposed by the condition that the partial wave is invariant 
under holomorphic and antiholomorphic conformal transformations of the coordinates z and z Jl2| . In terms 
of the wave functions x which depend only on parts of the coordinates, this condition leads to the following 
transformation properties 

x({z i: Zi};z ,z ) -> x({4> 4j z'o) = (cz Q + d) 2h (cz + d) 2h x({zi, z 2 }; z , z a ) (5) 

under conformal transformations z — > z' = az t b , an d z — > z' = aZi+b with ad — be = ad — be = 1. The possible 
values of the conformal weights h and h can be parameterized as follows [[l2| |J [| 

1 + m . j 1 - m 
h = — iv, h= — iv (6) 
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for arbitrary real v and integer m. 

It was suggested that these models can be exactly solved by means of the Bethe ansatz. Among other 
things, this would imply that the model contains a family of "hidden" integrals of motion. Indeed, at least one 
nontrivial integral of motion A was found It commutes with the hamiltonian and is given by 

A = z 12 z 23 ■■ -z nl Pi P 2 ---P n , [H n ,A} = (7) 

together with the analogous operator for the antiholomorphic system. It was noticed jjj], that this operator 
belongs to the family of n mutually commuting operators generated by a monodromy matrix which satisfies the 
Yang-Baxter equation. However it was not clear that all these operators commute with the hamiltonian and 
that they are the integrals of motion of the model. In the present paper we show that this is indeed the case. 
We prove that the model ([!]) is completely integrable and discuss the possibility of solving it by means of the 
Bethe ansatz. 



2. High energy multicolor QCD as XXX Heisenberg mag- 
net 

Let us show that the hamiltonians (||) and (Q) coincide with that of the XXX Heisenberg magnet for spin s = 0. 
To give meaning to the limit s = we recall that the quantum inverse scattering method allows one Jl3| to 
construct a family of exactly solvable XXX models for arbitrary complex values of the spin s. For arbitrary 
complex s the spin operators S k (a = 1, 2, 3) in all sites k = 1, . . . , n of the lattice can be realized as follows 

S+ = z 2 k d k - 2sz k , S k ~ = -d k , Sl = z k d k -s (8) 

The definition of the model is based on the existence of a fundamental matrix RfafaiX) which obeys the Yang- 
Baxter equation and acts in the auxiliary space Vf 1 <S> V/ 2 = h ® h with the space h having the dimension of the 



local quantum space in each site. For arbitrary complex spin s it is given by the following expression 14 



R hh (A) = f(s, A) r{lX _j niX + J+l) (9) 
where the operator J is defined in the space Vf 1 ® Vf 2 as a solution of the operator equation 

J(J + 1) = 25® S + 2s(s + 1) 

Notice that the definition of the R— matrix contains ambiguity in multiplication by an arbitrary c-number valued 
function /(s,A) which is fixed by imposing additional normalization conditions on Rf 1 f 2 (X). For the special 
value A = of the spectral parameter the operator Rf 1 f 2 (0) is proportional to the permutation operator Pf t f 2 on 
Yh. ® Vfa anc ^ both operators are equal to each other for finite dimensional representations of the spin operators 
p4[ provided that f(s, 0) = 1. The hamiltonian of the XXX model of spin s is given by the general expression 
(0) with the two-particle hamiltonian defined as 



#12 = \-^ lo & R hh( X ) 



(10) 



Let us consider the expression for the hamiltonian H in the limit s — 0. Using the expressions (Q) and (|T(J) we 
find that one of the T— functions in the numerator of (^) diverges as A, s — * 0. However we use the possibility 
to choose /(s, A) = T(iX + 2s + l)/T(iX — s) to avoid this problem and obtain the two-particle hamiltonian as 

H ik = -il>{-J ik ) - ^(J lfc + 1) + 2^(1) (11) 

Here ip(x) — d\ogT(x)/dx and the operator Ji k is one of the solutions of the equation 

Jik(Jik + 1) = 2Si <S>S k = -( Zi - z k ) 2 did k 



where we substituted the explicit form (|8|) of the spin operators for s = 0. Comparing the last relation with 

-2 



we notice that the operator Ji k {Ji k + 1) coincides with the definition of the operator L 2 - k . Substituting 
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Lj k = Jjk{Jjk + 1) into (j^j one can perform the summation over I to get an expression for the hamiltonian 
which is identical to that in ([ll])! Thus, we identify the holomorphic and antiholomorphic QCD hamiltonians 
([l]) as hamiltonians of a one-dimensional XXX Heisenberg model with spin s = 0. In what follows, we study 
only the holomorphic sector of the model. The generalization to the antiholomorphic sector is straightforward. 

The identification of the model as XXX magnet means that the system (Q), which describes high-energy 
asymptotics in multi-color QCD, is exactly solvable. To find the family of local integrals of motion of the model 
we follow the standard procedure. To any site k we assign auxiliary and fundamental Lax operators |l4| 

ifc,a(A) = XI„ <g> I a + iS k $ u a , L kJ {\) = R kJ (X) (12) 

where a are the Pauli matrices. The operators L kM and L k j act locally in the space ® <D 2 and hk <g) h, 
respectively, where hfc is the quantum space in the k— th site and the dimensions of h and hfc coincide. Taking 
the ordered product of the Lax operators along the lattice we define the auxiliary monodromy matrix as a 
matrix m <L 

T a (X) = L„ ia (A)L„_i !0 (A) . . . £*i, a (A) = ( c\\) -D(A) 1 

and analogously the fundamental monodromy matrix T/(A). Taking the trace of the monodromy matrices over 
the auxiliary space we get two operators, the auxiliary and fundamental transfer matrices, 

A(A) = tr„ T„(A), r(A) =tr/ T/(A) (14) 

which act in the full quantum space of the model and commute with each other for different values of the 
spectral parameters [|| 

[r(A), A(m)] = [r(A), r(/*)] = [A (A), A(/x)] = (15) 

Differentiating the both sides of this relation with respect to the spectral parameters A and fj, and putting 
A = /J, = Q we will get a family of mutually commuting conservation laws of the model. Moreover, the 
fundamental transfer matrix t(A) contains the local integrals of motion including the hamiltonian of the model, 
while the operator A(A) allows one to construct their eigenstates by means of the Bethe ansatz. The explicit 
form of the local integrals of motion is given by |14] 



1 d k 



— logtr.TMA) 



1,2,... (16) 



where the k— th operator describes the interaction between k + 1 nearest neighbors on the lattice. In particular, 
the operator I\ = H = £\ =1 H k .k+i coincides with the hamiltonian of the model. 

To find the explicit form of the operator A(A) we substitute the spin operators (ph for s = into the standard 
form of the auxiliary Lax operator ( p^ ) and represent the result as a matrix in C 

L k,7 0) W = XI + iv k ®u k d k (17) 

where d k = -t^- and v k = (g h J an d &k — ( z k, — 1) are vectors in the auxiliary space C 2 . After substitution of 
( |l7| ) into ( |l3| ) we use the relation u k Vj = (z k — Zj) to get 

A S=0 (A) = 2A n - S(S + 1)A"- 2 + Q 3 A"- 3 + . . . + Q n (18) 

where the operator S is defined as 

S(S + 1) = - J] z 2 kj d k dj = SS = S 3 S 3 + \{S + S- + S~S + ) 
k>j 

and S is the total spin of lattice. The operators Q k for k = 3, . . . , n are given by 

Qk ^ ^ 1 * Zi 1 i 2 Zi 2 i 3 . . . Zi k i 1 di 1 3i 2 . . . &i k 

n>ii>i2> ■ ■ ->ifc>l 
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and we recognize the operator Q n as the operator A defined in (|7]) . Notice, that in the expansion of A s= q (A) 
the term with A™" 1 is absent due to the orthogonality of vectors Vk and Uk- It follows immediately from ( |l5| ) 
that the integrals of motions Ik and Qj are mutually commuting operators 

[h,Q j ] = [Qk,Q j } = [h,I j ] = o 

with Ii = H being the holomorphic hamiltonian. 

To diagonalize these operators one may try to apply the algebraic Bethe ansatz However, the attempt 
to use the Bethe ansatz for s — fails from the very beginning. The algebraic Bethe ansatz is based on the 
existence of the highest weights in each cite of the lattice. Trying to define them, we find that for s = the 
equations S^\u>k) = and Sjl\uj) — s\ui) have a trivial solution \ui n ) = const which is annihilated by the spin 
operator S. The reason why it happens is the following. As we will show in sect. 4 the spin operators form for 
s = the representation t°' 2 (g) of the principal series of the noncompact group SL(2,€). This is a unitary 
infinite dimensional representation Jl5| which does not contain a highest weight, in contrast with unitary finite 
dimensional representations of SU(2). Thus, to diagonalize the hamiltonian of the model one has to generalize 
the algebraic Bethe ansatz to the case of noncompact groups. Nevertheless, there is a trick which allows one to 
apply the algebraic Bethe ansatz to find special eigenstates of the model. 



3. Algebraic Bethe ansatz for spin s = — 1 

There is one-to-one correspondence between the XXX models of spin s = and s = — 1 based on the following 
relation between the Lax operators (|l2|) in both models 

-(«=■ 

J k,a 



4 S : _1) (A) = XI + id k vk®u k = (L { k S - 0) (X)) T = P k L{ S - 0) (X) 



where operator T denotes similarity transformation. The fundamental Lax operators ( p^ ) satisfy the same 
relation. From this property and the definition (|l^) we find the relation between the monodromy matrices of 
the models 

Ti s= - 1] (X) = (Ti s =°)(A)) T = P 1 P 2 ...P n Ti s =°\X) (P 1 P 2 ...P n )~ 1 

and analogously for T^ s ~ ^(A). It is obvious from ( |l4| ) that the same relations are valid between the transfer 
matrices but in this case one can use the fact that they commute with the operator Q n = A defined in (Q) to 
write 

A s= _i(A) = Af =0 (A) = (zi 2 z 23 ■ ■ ■ z nl )~ x A S=0 (A) z 12 z 23 ...z nl (19) 

and analogously for r s= _i(A). Thus, the transfer matrices in the XXX models of spin s = and s = — 1 have 
the same eigenvalues and their eigenstates are related as follows 

\<p(zi', z Q )) = z 12 z 23 . . . Zni\<p(zi; z )) (20) 

As follows from (|l0|) the same relation (jl^) holds between the hamiltonians of both models. 

The reason why we included the XXX model of spin s = — 1 into consideration is that one is able to apply 
the algebraic Bethe ansatz for spin s = — 1.Q In contrast with spin s = 0, it is possible to find a nontrivial 
highest weight for spin s = — 1. Indeed, by substituting the explicit form of the spin operator (||) for s = — 1 
into the equations St\uJk) = and Sf.\uJk) = —\^k) w e find that it has the nontrivial solution \u>k) = \j z\ which 
allows one to construct the pseudovacuum state as 

1^ = 33^ 



=1*2 



To define the Bethe states one has to use the B operator for spin s — — 1 defined in (|T^). Then the Bethe states 
are given by 

= B(X 1 )B(X 2 ) . . . fl(A ; ) } (21) 



1 However, as is shown in sect. 4, the algebraic Bethe ansatz gives us only special set of the eigenstates 
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where the parameters {A} = (Ai, . . . , A;) are solutions of the Bethe equation for s = — 1 



= n w 

and the operator B is defined through ( |l3| ) as an element of the auxiliary monodromy matrix corresponding 
to spin s = — 1. The Bethe states ( pl| ) depend on a positive integer number I which is fixed by the additional 
constraint (§). Indeed, using ( p0| ) it can be easily shown that the transformation properties (||) of an eigenstate 
of the XXX hamiltonian for s — are equivalent to the following conditions for the eigenstates of the XXX 
hamiltonian for spin s = — 1 

S 3 |fe0)) = 0)), S+|£(z i; 0)> = 0, \<p{ Zi ; z )> = \<p( Zi - z : 0)) (23) 

where S is the total spin of the lattice for s = — 1 . Comparing these constraints with the properties of the Bethe 
states, S 3 \0) = — (n + t)\<p) and S + \<p) = 0, we find that 

n + l^h, I = 0, 1, .. . 

where h was defined in (|J) and n is the number of sites, or equivalently, the number of reggeized gluons in the 
t— channel. This relation means that the algebraic Bethe ansatz allows us to construct a very special class of 
eigenstates corresponding to positive integer values of the conformal weight h of the n— gluon wave function. 
To find the remaining states for arbitrary complex h a generalization of the Bethe ansatz is required. 

Using the relations ( |2l| ) and ( pp| ) we finally get the expression for the Bethe states which diagonalize the 
original hamiltonian as follows 

\(fii(z z ; 0; {A})} = 212^23 ■ ■ • z n iB(\i)B(\ 2 ) . . . B(X t ) „ „ 1 ? 



■n*2 • • - A n 

The integrals of motion for the XXX model with spin s = — 1 are related to those for XXX model with spin 
s = by equations identical to equation (|l^) between the transfer matrices. That is why they have the same 
eigenvalues in both models. The explicit expressions for the eigenvalues of integrals of motions for arbitrary 
spin s have been found in algebraic Bethe ansatz jl4| and we use these expressions for s = —I to get 

^E-^a^- (24) 

3=1 J J 

Note that this expression was found with the function /(s, A) that defines the normalization of the R— matrix 
in (||) chosen to be equal to one. It follows from ([If]), that for arbitrary /(s,A) the expression ( p4| ) gets an 
additional trivial contribution — in-jrzf(s, A)|a=o- Given a solution to the Bethe equation ( p2|) this relation 
yields the spectrum of local integrals of motion in the original model with hamiltonian H . In particular, the 
eigenvalues of the hamiltonian are equal to 

1 1 



£n = J 1 (A) = -2^_ I 



k=l 



where {A} obey the Bethe equation ( |22| ) for a fixed number of reggeized gluons n and s = —I. Thus, for special 
values of the conformal weights h the algebraic Bethe ansatz constructed in this subsection yields the wave 
function of the model ([!]) and the corresponding eigenvalue of the QCD hamiltonian. 

4. Eigenstates as highest weights of representations of 

SX(2,(D) 

So far we have considered only the properties of the holomorphic hamiltonian in (|l]) . If we know its eigenstates 
we can generalize them for the antiholomorphic hamiltonian in (Q) and construct the solution of the system ([j]) 
as the product of holomorphic and antiholomorphic eigenstates 

X({zi, Zi}; z , zq) = ip{zi\ z )<p(z~i] zq) 
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Let us combine the holomorphic and antiholomorphic hamiltonian to define the following operator acting on 
both z and z coordinates of particles: 

Ti = Bj % + H n 

where similarity transformation was defined in (|l9|). This hamiltonian has the following properties. As we 
have seen in sect. 3, the hamiltonian H n describes the XXX magnet with the generalized spin s = while the 
hamiltonian corresponds to the XXX magnet with spin s = — 1. Recall that the hamiltonians H n and H% 
are related to each other by an equation similar to ( [l9| ) and they have the same spectrum of eigenvalues and 
the corresponding eigenstates satisfy relation ( pp| ) . Hence, for fixed e n and e n in ([!]) the hamiltonian Ti has the 
eigenvalue e n + e n and the corresponding eigenstate is 

X({zi, Zi}; z , z ) = ip(zi\ z )ip(zi; z ) = (z 12 z 23 ■ ■ ■ z nl )~ 1 x({z l , £;}; z , z ) (25) 

Thus, there is one-to-one correspondence between this function and the eigenstate X- We notice that the 
hamiltonian Ti is selfadjoint and its eigestates x are orthogonal to each other. The reason why we consider the 
functions x({zj, z i\\ z o 5 z o) is that these functions are the highest weight of the irreducible infinite dimensional 
unitary representation of the principal series of the SL(2,C) group. 

To show this we recall that the unitary representation of the principal series, t p,m (g), is realized on the space 
of square integrable functions f(z,z) according to 

t p > m (g)f(z,z) = (cz + d)- k l 2+i ^-Hcz + d) k ' 2+i Pl 2 - l f( ^1±1 

\cz + a cz + d 

where k is an integer, p is a real number and the scalar product in this space is given by 

Lfi\h) = \ dzdzfi(z,z)f%(z,z) 



In this representation the group generators can be realized as holomorphic spin operators (^) with s = — 1/2 — 
k/4 + ip/4 and antiholomorphic spin operators with s = — 1/2 + k/4 + ip/2. 

For the one-dimensional model with the hamiltonian Ti we have in each site k of the lattice holomorphic 
spin operators and S| defined by (§) for s = — 1 and antiholomorphic spin operators S^ 1 and 5| given by 
expressions d§) for s = with the replacement of z by z. These six operators commute with the hamiltonian Ti. 
(and not with the QCD hamiltonian H n + H n l) and form the representation t ' 2 of the principal series of the 
SX(2,C) group which we identify as a local quantum space in the k— th site. The full spin operators of the lattice 
S a = J2k=i an d "5" — J2k=i ^k act 011 t ne tensor product of these representations which is reducible with 
respect to SX(2,C) and which can be decomposed into the direct sum of the irreducible unitary representations 
of the principal series as follows. It is well known Esl that the product of two infinite dimensional unitary 
representations of the principal series, t pi ' 1 and t P2 ' 2m2 : where mi and m2 are integer and pi and p 2 are real, 
is decomposed into the direct sum of the irreducible representations of the same principal series t p3 ' 2 with 
m,3 being integer. Applying this result to the tensor product of n copies of i ' 2 corresponding to all site of the 
lattice we get 

n P' m 

that is the total quantum space of the lattice is decomposed into the sum of infinite dimensional unitary 
representations of the principal series of SX(2,C) with real p and integer m. In principle, one may try to evaluate 
the multiplicities and the Clebsch-Gordan coefficients for each t p - 2rn in this decomposition using the standard 
formulas ]l5[ . Thus, the eigenstates x{{ z ii zq, zq) of the hamiltonian Ti belong to the representations of 
the principal series t p,2m . However, there are additional constraints on the functions x({^i? -^}; zq, Zq) imposed 
by the transformation properties (||) , or equivalently by the conditions ( p3| ) for holomorphic states and similar 
relations for antiholomorphic states. Using the relation (^) we represent these conditions in the following form 

S 3 \x) = -h\ X ), S+\x) = 0, S 3 \x) = -h\x), S+\ X )=0 

where x = X.{{ z ii z i}i z o = 0, zq — 0) and the conformal weights h and h were defined in ([|). Using the explicit 
form of h and h we conclude that the eigenstates of the hamiltonian TC which obey these additional conditions 
are the highest weights of the irreducible representation t il/ ' 2m of the principal series of the SL(2,<C). Thus, 
for arbitrary complex conformal weights h and h in (^) there might exist the eigenstate of the model. On the 
other hand, as was shown in sect. 3, the algebraic Bethe ansatz is applicable only for integer h and it cannot be 
generalized even for noninteger conformal weights. 



7 



5. Generalized Bethe ansatz 



To find the eigenstates corresponding to arbitrary complex values of the conformal weights one could use the 
method of the Q— operator [jl6| proposed by Baxter. Using the relations between the XXX models of spins 
s = and s = — 1 we restrict ourselves only by consideration of the model with s = — 1. In this case, there 
exists an operator Q(X) which acts on the full quantum space of the lattice and obeys the Baxter equation 

A(A)Q(A) = (A - i) n Q{\ - i) + (A + i)"Q(A + i) , [Q(A), A(/*)] = [Q(A), Q(/i)] = (26) 

Thus, the operators A and Q can be diagonalized simultaneously. Their corresponding eigenvalues are c— number 
valued functions of the spectral parameter A which obey the same equation ( |26| ) . That is why we are using below 
the same notations A(A) and Q(X) for the eigenvalues of these operators. Solving this functional equation and 
using the analyticity properties of A(A) as a function of A, one is able to find the eigenvalues of the auxiliary 
transfer matrix A (A) and the operator Q(X). The standard way of solving the Baxter equation consists of 
choosing eigenvalues Q(X) to depend polynomially on the spectral parameter A 

/ 

Q(A) = const. JJ(A-A fe ) (27) 
fc=i 

where {A} = {Ai,...,A;} are arbitrary complex numbers and / is positive integer. The parameters {A} are 
found by imposing analyticity conditions for the function A(A) at A = A&. Namely, substituting A = A^ into 
the both sides of the Baxter equation we get an equation for the parameters {A} which coincides with the 
Bethe equation (^2|). Thus, the algebraic Bethe ansatz corresponds to the special case of the eigenvalues of 
the Q— operator as finite polynomials. To go beyond the algebraic Bethe ansatz one has to look for solutions 
of the Baxter equation among the functions more general than polynomials. The simplest possibility which is 
realized in the Toda model is to search for the solution among the entire functions having infinite number of 
zeros, I — » oo. 

Using the definition ( |l8| ) we find that the function A(A) is a polynomial of power n in A with the coefficients 
given by the operators S(S+ 1) and Qk- This suggests the following way to solve the Baxter equation. First, one 
substitutes the general form of the function A(A) into the l.h.s. of the Baxter equation, replacing all operators 
by their unknown eigenvalues, S and Qk (k = 3, . . . ,rt). Second, one solves the functional equation ( p6| ) and 
finds the solution for Q(X) which depends on these quantum numbers. Third, imposing the condition that Q(X) 
is an entire function of A one determines the possible values of 5* and Qk and substitutes them into ( |l8| ) to find 
the eigenvalue of the auxiliary transfer matrix. 

Thus, for any fixed allowed set of the quantum numbers {S, Qk} we will find the corresponding function Q(X) 
satisfying the Baxter equation. It is important to notice that among these functions there are very special ones 
( pTj ) corresponding to the Bethe ansatz solution of the XXX model of spin s = — 1. Moreover, in this particular 
case we know the explicit form of the eigenstates (|2l]) and eigenvalues ( pi| ) of the local integrals of motion. All 
these expressions, including the Q— function (p7j), explicitly depend on the parameters {Afe} which are solutions 
of the Bethe equatio n (j22"|) . The natural question appears: is it possible to express both the eigenstates ( f2lj ) 
and the eigenvalues (|24[)in terms of Q(X) in such a way, that all their dependence on {Xk} will be contained 
inside the Q— function? In this case, one will get a unique possibility to perform an "analytical continuation" 
of the results found within the framework of the algebraic Bethe ansatz by replacing the special values (|27|) of 
the Q— function by a general solution of the Baxter equation. 

Let us express the Bethe states (|2l| ) in terms of the Q— function defined in (p7j). Substituting the explicit 
form of the auxiliary Lax operator into the definitions (^) one finds that for s = — 1 the operator B(X) is a 
polynomial of power n — 1 in A which can be represented in the following form 

B(X) = iS~(X - xi)(A - x 2 ) . • • (A - x n -x) (28) 

where x\, . . . , x n -\ are operator zeros of B(X). The precise meaning of Xk as roots of the operator polynomial 
B(X) was given in the framework of the functional Bethe ansatz |l7], These operators have the following 
properties 

[S~,x k ] = [S 3 ,x k ] = [xk,Xj] = 



which allows us not to worry about their ordering in feq). Now we may rewrite the Bethe states (x2W) in terms 
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of the operators Xk 

l n—1 n — 1 I 

\h) = n is ~ - x ^ = ( is ~y n ri( A i - 
j=i k=i k=i j=i 

where Xj satisfies the Bethe equation (^4) for s = — 1 and fixed I and n. Now we immediately recognize that 
the product over j coincides with the expression (^7j) for the function Q(A) defined for A = Xk- Notice that 
the Q— operator doesn't commute with the operators Xk and the notation Q(xk) means the eigenvalue of the 
Q— operator defined to be evaluated for operator value of the spectral parameter. Thus, the Bethe states can 
be expressed in terms of the eigenvalues of the Q— operator as follows 

|0) = (iS-) h - n Q(x 1 )Q(x 2 ) . . . Q( Xn - 1 )- s -^ (29) 

z l z 2 ■ ■ ■ Z n 

where h is the conformal weight defined in (Q). The remarkable property of this expression is that it does not 
depend explicitly on the parameters A& . The operators S~ and Xk are determined by the properties of the model 
and are the same for all eigenstates. Different eigenstates \4>{) are parameterized by different solutions <5(A) of 
the Baxter equation. After being written in this form, the relations (^9|) admit a natural analytic continuation 
to arbitrary complex values of S, or equivalently h. Moreover, these states satisfy the additional conditions ( p3| ) 
which follow from the analytical continuation of the analogous relations for the Bethe states. 

Once we know the eigenvalues of the Q— operator the eigenvalues of the auxiliary transfer matrix can be easily 
found from (|2(]). To find the eigenvalues of the local integrals of motion we use their expressions found by the 
algebraic Bethe ansatz method, rewrite them in terms of the Q— functions and perform analytical continuation. 
Using equation (|24|) we represent the eigenvalues of the operator for s — — 1 in the following form 

^--^^^—--d^^ll Xj -i + X 

3 = 1 3 J 0=1 J 



(30) 



and in terms of Q— function fl27| ) this expression looks like 

1 d k Q(-A-i) 
7 *= ^-7u lo S 



idX k °Q(-X + i) 



Although this relation was found as a relation between eigenvalues of the operators Ik and Q it can be extended 
to be an operator relation because both operators have the same eigenstates. 



6. Solution of the Baxter equation 

To find the eigenstates and eigenvalues of the integrals of motion one has to solve the Baxter equation for 
s = —1 and substitute the solution for the function Q(X) into the relations (|2^) and (|30|). Using the general 
form of the function A(A) in terms of quantum numbers {S, Qk} we get the following functional equation for 
the function Q(X) 

(2A n - h(h - l)X n -' 2 + Q3A"- 3 + . . . + Q n )Q(X) = (A + i) n Q{\ + i) + (A - i) n Q(X - i) (31) 

where we used the relations (p3| ) and (|l^) to identify the eigenvalue of the operator S with the conformal weight 
h of the eigenstates defined in (^[). In this equation h and Qk are independent parameters whose values are 
restricted by analyticity properties of the solution Q(A). 

The Baxter equation ( |3l| ) depends on the set of quantum numbers {h, Qk} and obeys the following property. 
It is invariant under the replacement 

h -> 1 - h (32) 

Hence, if the function Q(X; h, {Qk}) is a solution of the Baxter equation then so is Q(X; 1 — h, {Qk})- Combined 
with (|29l), this property allows us to relate the eigenstates of the model with the conformal weights h and 1 — h. 
Thus, trying to solve the Baxter equation ( ^l|) we use the symmetry (|3^) to restrict the possible values of the 
conformal weights to the fundamental region 

Re h > - (33) 
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or in terms of the parameters m denned in (|6j), m > 0. It is interesting to note that the Bethe ansatz solution 
( pTj ) for 0(A) corresponds to integer positive h and the symmetry ( |32| ) allows us to generalize the results of 
the sect. 3 to arbitrary integer h by putting I = \h — ^| — 



i. This is the simplest example of the analytical 



continuation of results obtained by means of the Bethe ansatz. 
In the limit of large A the Baxter equation can be rewritten as 



-h(h - l)X n - 2 Q(X) = (A + i) n Q{\ + i) + (A - i) n Q(X - i) - 2\ n Q{\) ~ ^_(A"Q(A)) 

where we neglected nonleading terms in both sides of the equation. Solving the resulting second order differ- 
ential equation we get the general large— A asymptotic behavior of the solutions of the Baxter equation in the 
fundamental region ( |33| ) 

A' (34) 



0(A) A A' 1 - 



Using the asymptotics ([mJ) we recognize the special role of integer positive values of the conformal weights h. 
The analyticity of 0(A) as a function of A together with (34) implies that for h — n G 2 + the solution of the 
Baxter equation is polynomial of order h — n in A which is in accordance with ( p7j ) . 

To solve the Baxter equation (J3l|) for arbitrary A we perform the Mellin transformation of the function Q(A) 



0(A) = 



duj 



OH 



(35) 



After substitution of (£55|) into (|3l]), the original finite differential equation for 0(A) is replaced by an ordinary 
ri— th order differential equation for Q{lo): 



duj 



k=0 



dio 



Q{u) = 



where Q2 = —h(h — 1). A simple analysis shows that the solution of this equation has singularities at u> = 1. 
Indeed, in the large A limit the function (|35| ) gets its leading contribution from the integration at the vicinity 
of the point ui — 1 and the asymptotics (|34|) imply that 



OH "^(i-wy 



l-h+n 



where h takes values in the fundamental region (^). Thus, for the expression ( p5| ) to be well defined one has to 
fix the prescription for the integration of the lj = 1 singularity. Since the integration in (|3^) is performed along 
the real positive axis in the complex to— plane one can deform the integration path at the point u) = 1 to encircle 
the singularity either in lower or upper half plane. Moreover, rotating both integration paths as lo — > e lw uj and 
lu — > c~ l7r w, respectively, we avoid the singular point lo — 1 and get the following expression for the Q— function 



/>oo 

Q(A) = (C+e wX + C_e" 7rA ) / duj 

Jo 



Q(-lj) 



where C+ and C_ are arbitrary constants. However the value of the constants is fixed by the condition that for 
arbitrary complex A the integral should be convergent at to — and w = 00 



where the numerical value can be arbitrary. Changing the integration variable z = 1/(1 — lu) we obtain 



(36) 



Q(A) = sinh(TrA) / dz (1 - z _iA_1 Q(z) 



(37) 



where the function Q(z) obeys the following equation 

n-2 



(1-*)^) +z{l-z)Y,Qn- k (-lz(l-z)- 



k=0 



dz 



Q{z) = Q 
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Let us consider the solutions of this equation in the special case n = 2 corresponding to the compound state of 
two reggeized gluons - the BFKL pomeron. 

For n = 2 the equation for the function Q(z) has the following form 

z(l - z)Q"(z) + (1 - 2z)Q'{z) + h(h - l)Q(z) = 

where prime denotes differentiation with respect to z. The equation is invariant under the replacement z — > 1 — z 
and its general solution is well known [fL9| to be a linear combination of the hypergeometric functions 

Q(z) = aF(h, 1 - h; 1; z) - c 2 F(h, l-h;l;l-z) 

where c\ and C2 are arbitrary constants. After substitution of this expression into ( |37| ) we find the Q— function 

as 

Q(X) = ciQo(A) + c 2 Q (-A) (38) 

where Qo(X) denotes the function 

Qo(A) = sinh(TrA) / dz (1 - z)^' 1 z~ lX ~ x F(h, 1 - h; 1; z) 



We use the integral representation for the hypergeometric function ]i9| | in order to express the solution in the 
following form 

Q (A) = sinh(TrA) Bmfyh) f dz (1 - z )'A-i z -iA-i f dt t h-l^ _ t yh^ _ ^h-i ^ 

Jo Jo 

This expression for Qo(X) can be represented as a double Pochhammer contour integral 

Q (A) = dzz-^iz-l)^- 1 j dtt h - 1 (t-l)- h (l-tz) h - 1 

where the contour P incloses the singular points and 1 in complex z and t planes. Taking (|3^) and expanding 
the last factor in the integrand in powers of zt we perform the integration over z and get the following factor 
r(zA)r(l — iX) = — gfrj^yjg which generates the singularities of the function Q(X) at A = ilL. They originate 
from the "dangerous" points z = and z = 1 in the integral. To preserve the analyticity of the Q— function we 
use the possibility to choose arbitrary values of the constants C+ and C_ to put them equal to ([36]) . This leads 
to the appearance of the factor sinh(7rA) in (^) which compensates the singularities of the integral. Finally, we 
get the representation for Qq as an infinite sum 

where is the generalized hypergeometric function Jl9| . This expression has the following properties. 

In the k— th term of the sum the A— dependence comes from the ratio T(k — iA)/r(l — iX) which is a 
polynomial of order k in A. Hence, being expanded in powers of A the function Qo{X) turns out to be an infinite 
series. This seems to be in contradiction with the fact, that in the special case of positive integer values of 
I = h — 2 the Bcthc ansatz gives us an expression for the Q— function which is a polynomial of power h — 2 in 
A. However, notice that in the representation j40| ) the factor 1/T(h — k) ensures the truncation of the sum after 
the k = (h — 1)— th term for h > 2. Moreover, the expression ( [40| ) explicitly obey the symmetry ( |32| ) which 
implies that the series terminates also for negative integer h. 

For arbitrary noninteger h the expression (^) is an infinite series in A. Using the properties p9[ of the 
function 3F2 one finds that this series converges only if Im A < 0. Then, the general form of the solution ( p8| ) 
implies that the solution of the Baxter equation Q(A) is given by the function Qo(X) in the lower half plane in 
A and by Qo{— X) in the upper half plane. It is interesting to note that for integer h the function satisfies the 
relation Qo{— A) = (— l) h Qo(X) which implies that both terms in ( |38|) are equivalent. 

After substitution of ( fiol) into ([38]) we get an expression for Q(X) which is an entire function of A having an 
infinite number of zeros. Only in the case of positive integer h the number of zeros is finite. It can be checked 
that the roots {Afc} do obey the Bethe equation for n = 2. This is in accordance with general conditions on the 
Q— function discussed in sect. 5. 
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7. Conclusions 



In this paper we found that the one-dimensional lattice model proposed by Lipatov to describe high-energy 
scattering of hadrons in QCD is completely integrable. Applying the quantum inverse scattering method we 
identified the Lipatov model as the generalized one-dimensional XXX chain of spin s = 0. We found the 
family of local integrals of motion and developed the generalized Bethe ansatz for their diagonalization. The 
corresponding eigenstates and eigenvalues, (^9|) and (|30|), are expressed in terms of the Q— function which 
satisfies the Baxter equation (^6|) for spin s = — 1. Solving this equation in the special case of the lattice with 
n = 2 sites we found the expressions ( |38| ) and ( f40| ) for the Q— function. It turns out that after the substitution 
of this function into ( p9| ) and ( |30| ) we will get the expressions for the wave function \ of the compound state of 
n = 2 reggeized gluons and the corresponding eigenvalue of the QCD hamiltonian which are identical to that 
|jl2[ for the BFKL pomeron. The details of the calculations and the proof of the statements made in this paper 
will be published elsewhere. 

We would like to thank V.Korepin and L.Takhtajan for helpful conversations. L.D.F. is grateful to C.N.Yang 
for the hospitality at the Institute for Theoretical Physics at Stony Brook which in particular made our collab- 
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for their support. The work of G.P.K. is supported in part by the National Science Foundation under grant 
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